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ABSTRACT 

We present numerical simulations that include 1-D Eulerian multi-group radiation- 
hydrodynamics, 1-D non-Local-Thermodynamic-Equilibrium (non-LTE) radiative transfer, 
and 2-D polarised radiative transfer for super-luminous interacting supernovae (SNe). Our 
reference model is a ^IOMq inner shell with 10®^ erg ramming into a ~3M0 cold outer 
shell (the circumstellar-medium, or CSM) that extends from 10^^ to 2x10^® cm and moves at 
100 km s“^. We discuss the light curve evolution, which cannot be captured adequately with a 
grey approach. In these interactions, the shock-crossing time through the optically-thick CSM 
is much longer than the photon diffusion time. Radiation is thus continuously leaking from the 
shock through the CSM, in disagreement with the shell-shocked model that is often invoked. 
Our spectra redden with time, with a peak distribution in the near-UV during the first month 
gradually shifting to the optical range over the following year. Initially Balmer lines exhibit 
a narrow line core and the broad line wings that are characteristic of electron scattering in 
the SNe Iln atmospheres (CSM). At later times they also exhibit a broad blue shifted com¬ 
ponent which arises from the cold dense shell. Our model results are broadly consistent with 
the bolometric light curve and spectral evolution observed for SN 2010jl. Invoking a prolate 
pole-to-equator density ratio in the CSM, we can also reproduce the ~2% continuum polarisa¬ 
tion, and line depolarisation, observed in SN 2010jl. By varying the inner shell kinetic energy 
and the mass and extent of the outer shell, a large range of peak luminosities and durations, 
broadly compatible with super-luminous SNe Iln like 2010jl or 2006gy, can be produced. 
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1 INTRODUCTION 


Since the or i ginal i 
IPopita et"^ 11984 


Niemela et al.l 

kO 

00 

Schlegell ll99C 

), the 


culiar Type II SNe has grown from a few to a few tens 
of events. It has also revealed an intriguing diversity, span¬ 
ning a range of luminosities, durati on, line profile morphol¬ 
ogy, with events like SNe 1988Z l Stathakjs&Sadlgll I99l|; 


2004h. 2006gv 

Smith et alj 12007), and the enigmatic SN2009ip 

I Mauerhan et al 

I 2 OI 3 I; Pastorello et al.l2013l:lMargutti et alj20l4. 


The origin of Type Iln SNe is associated with the interaction of 
fast material with previously -expelled slowly-expanding material 
JOrasberg & Nadezhinll98al . The resulting shock leads to deceler¬ 
ation of the inner shell, acceleration of the outer shell, and con¬ 
version of kinetic energy into internal and radiative energy. Be¬ 
cause (standard) non-interacting SNe radiate merely 1% of their 
total energy (the rest being kinetic), interacting SNe can easily 


reach extraordinary luminosities by tapping this abundant reser¬ 
voir of kinetic energy. For example, an interaction extracting 30% 
of a 10®^ erg ejecta could produce a SN Iln with a luminosity of 
10^° Lq for 3 months. 


Super luminous interacting SNe require both a high kinetic 
energy in the i nner shell and a larg e mass reservoir in the outer 
shell (see, e.g.. lMoriva et ^ 201_^). T wo go od examples of such 
SNe Iln are SN 2006gy_j Smith et ^ l2007h and more recently 
SN2010jl jStoll et aljEoulT The latter has extensiv e photomet¬ 
ric and spectroscopic observations i Stone^i]|201 ll : IZhang et al.l 
I2OI2I : lOfek et ai] I2OI4 ; Frans son et alj 201 4 and also has spec- 


tropolarimetric data l lPatat et al.ll201 iLIWilliams et al.ll20l4) . How 
stars produce the outer shell is unclear - it may be mate¬ 
rial ejected by the p air-production instability jBarkat et al .11 19671 : 
IWooslev et al.l |200% or it may be material ejected via a super- 
Eddington wind jOwocki et alj2004h . Here, we simply consider the 
interaction between a fast moving inner shell and a slowly-moving 
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massive outer shell, and leave aside speculations concerning their 
origin. 

Numerical simulations of interacting SNe are generally 
undertaken with radiation hydrodynam ics jChugai et al.l l2 004l : 

IWooslev et al.ll200^ : iMoriva et al.ll2013^ : IWhalen et al.ll2013l) , hy¬ 

drodynamics combined with a parame terised cooling function 
jvan Marie et alJl201(ll: Ehen et alJl2014h. or pure hydrod ynamics 
(hy assuming adiabaticity; see, e.g., Blondin et all 1 199m . When 
studying super-luminous SNe Iln, radiation hydrodynamics with 
a non-grey (multi-group) treatment is a must because the radia¬ 
tion contribution is a sizeable fraction of the pressure/energy in the 
problem, the characteristic temperatures of the radiation and of the 
gas may differ strongly, and the emerging radiation provides essen¬ 
tial signatures characterising the event. 

The radiative transfer modelling for SNe lln comes at present 
in two flavours. It is either based on the results of radiation hy¬ 
drodynamics simiflations but limited to a single line, e.g., Ha 
dChugai et al.ll2004h . or based on an ad-hoc atmospheric structure 
and the assumption of steady -state radiation thro ugh a diffusive op¬ 
tically thick inner boundary jPessart et alj200^ . 

In this work, we discuss numerical simulations for super- 
luminous SNe of type Iln, including studies with radiation hydro¬ 
dynamics, radiative transfer, and polarised radiative transfer codes. 
In Section]^ we first present the multi-group 1-D Eulerian radia¬ 
tion hydrodynamics simulations we perform with HERACLES, in¬ 
cluding the numerical approach, the standard interaction configu¬ 
ration we adopt as initial conditions, and the results for both the 
dynamics and the radiation. We then present in Section[3]non-LTE 
radiative-transfer simulations with CMEGEN, which we compute 
for a number of snapshots using the results from the radiation- 
hydrodynamical modelling. We discuss line profile broadening, and 
the relative roles of expansion and electron scattering. In Section|4] 
we present polarisation calculations based on these CMEGEN simu¬ 
lations but artificially distorted along an axis of symmetry, in order 
to produce prolate density configurations. We assess the level of 
asymmetry needed to explain the observed continuum polarisation 
(and line depolarisation) reported for SN2010jl. In Section |5] we 
discuss the dependency of our results on the interaction properties, 
varying the kinetic energy of the inner shell and the mass/extent of 
the circumstellar medium (CSM). Einally, in Section|^ we present 
our conclusions. 


2 RADIATION HYDRODYNAMICS WITH Heracles 
2.1 Numerical approach 


HERACLES is a Eulerian multi-dim ensional radiation- 
hydrodynamics code dOonzalez et al.l 120071), with t h e pos sibility 
for multi-group radiation transport l Vavtet et al.l 1201 ih . The 
hydrodynamics is treated using a standard second order Godunov 
scheme. For the radiation transfer, con duction, flux limited diffu- 
sion (ELD) and the Ml moment model l lDubroca & Feuge^l 19991) 
are implemented in HERACLES, the latter two with the possibility 
of using a multi-group approach. Apart from a comparison test that 
uses a grey treatment, the simulations presented in this paper are 
all done using the multi-group Ml moment model. As explained 
below, a multi-group treatment is necessary in order to take into 
account the large opacity variations with frequency and the strong 
imbalance between radiation and gas temperature. The Ml model 
has also some strong benefit compared to the more standard 
ELD approach. First of all the Ml model is exact in both the 


diffusion and the free-streaming limit and is rather accurate in 
between. It is therefore well suited to model regions with low or 
intermediate optical depth. Furthermore, the Ml method solves 
explicitly for both the energy and the flux through moments of 
the radiative transfer equations. The Ml model can distinguish 
between absorption/emission and scattering opacities. The effect 
of scattering in regions where the scattering opacity is dominant 
can then be taken into account. Finally, the radiativ e transfer is 
formu lated in the co-moving frame and as is shown in IVavtet et al.l 
( l201lh the Doppler and aberrations effects are properly accounted 
for, which is important when large velocity gradients are present. 

We use the multi-group solver because grey radiation trans¬ 
port, often used for simulations of non-interacting Type II SNe, is 
inadequate in SNe Iln. In non-interacting Type II SNe, the radiation 
and the gas have similar temperatures everywhere (even identical at 
large optical depths), so that adopting the Rosseland mean opacity, 
for example, is not a bad representation of the effective opacity of 
the material. In contrast, the radiation produced from the interac¬ 
tion in SNe Iln has a characteristic temperature that can vastly differ 
from the medium that it traverses, and the total continuum optical 
depth is not much in excess of about ten. In the next section, we il¬ 
lustrate the impact on the light curve produced by assuming a grey 
rather than an energy-dependent LTE opacity for the material. 

In this work, we adopt a uniform composition, with a H mass 
fraction of 0.633, He mass fraction of 0.36564, and an iron mass 
fraction of 0.00136 to reflect approximately the near-solar compo¬ 
sition of a blue or a red su pergiant (hereafter BSG/RSG) star (for 
SN 2010iL [^oll et alj201lh argue for a sub-solar metallicity, which 
may be important for understanding how these SNe Iln come about, 
but is irrelevant for the present radiation-hydrodynamics consider¬ 
ations). 

We supply the code with an opacity table for our adopted com¬ 
position. We compute our opacities as a function of density, tem¬ 
perature and energy group. Energy groups are positioned at strate¬ 
gic locations to capture the strong variation in absorptive opacity 
with wavelength. We use one group for the entire Lyman contin¬ 
uum (including the X-ray range), two groups for the Baimer con¬ 
tinuum, two for the Paschen continuum, and three groups for the 
Brackett continuum and beyondQ Our opacity code computes the 
LTE level populations and ionisation state of the gas, and then uses 
the atomic data available within CMEGEN to compute LTE opaci¬ 
ties and to make the opacity table. The domain covered is 10 ~^° to 
10“® g cm“^ in 200 bins, and from 2000 to 500,000 K in 500 bins 
— we use the nearest edge for requests outside the table bound¬ 
aries. Here, we consider the contribution from electron scattering 
and bound-free opacity for H I, He I- H, and Fe I-XII. A major lim¬ 
itation of our work is that we ignore line opacity and emissivity 
in our HERACLES simulations. Metal line blanketing (in particular 
from Fe) would enhance the opacity in the UV (inhibiting further 
the escape of UV photons). We also neglect cooling due to CNO 
elements, which requires a non-LTE study with specific allowance 
for enhanced photon escape due to the velocity field. Inclusion of 
these effects is left to a future study. 

In Fig. [T] we plot the photon mean free path and its inverse 


^ We tried using several energy groups to cover the X-ray and the EUV 
ranges better but this caused numerical problems with the radiative trans¬ 
fer. The solver struggled converging when an exceedingly small energy/flux 
was present in those energy groups, a situation encountered outside of the 
shocked region. In practice, this should not impact our results sizeably be¬ 
cause the optical depth to hard radiation (all Lyman photons) is huge at all 
times considered here. 
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Figure 1. Top row: Illustration of the variation of the photon mean free path (left) and inverse mean free path (right) at a mass density of 1.16x 10“^"^ gcm“®. 
The electron-scattering opacity (black) is the main contribution to the Rosseland-mean opacity (red; which accounts for bound-bound and bound-free pro¬ 
cesses), except in the regions of partial ionisation of H and He at 8000-20000 K and in the iron-opacity bump in the range 100000-200000 K. The precise 
temperature regions where these offsets occur depend on density. The most critical variation in the opacity is associated with the recombination of H and 
He, leading to a very small opacity at low temperature. Bottom row: Illustration of the wavelength dependence of the photon mean free path associated with 
absorptive processes (blue) at ~ 3000 K (left) and ~ 12000 K (right), corresponding to conditions in which hydrogen is neutral and ionised, respectively. The 
filled dots con'espond to the central wavelength of each energy group used for both the opacity table and the HERACLES simulations. 


(i.e., the opacity) against wavelength for different processes at dif¬ 
ferent densities and temperatures. In super-luminous SNe Iln the 
spectrum formation region is typically optically thick in the Ly¬ 
man continuum (and at X-ray wavelengths) and high-energy pho¬ 
tons will tend to thermalise, even at low electron-scattering opti¬ 
cal depth. Longward of the Lyman continuum, electron scattering 
opacity typically dominates over absorptive opacities and so low- 
energy photons will thermalise inefficiently at most times and lo¬ 
cations. 

The simulations in this work use an ideal gas equation of state 
with 7 = 5/30 We also adopt a mean atomic weight of 1.38 in all 
simulations. In practice, this value depends on the ionisation state 
of the gas, but in ways that need to be determined by using a general 
equation of state — this is left to a forthcoming study. 

Our simulations are 1-D and use a uniform radial grid with 
1600 points, covering from 0.05 to 2x 10^® cm. Degrading the res- 


^ We therefore ignore the variations in thermodynamic quantities that arise 
from variations in ionisation. 


olution by a factor of four produces the same overall properties 
(e.g., the bolometric maximum is changed by ~ 10%), although the 
shock is less resolved and more numerical diffusion occurs as the 
interaction region crosses grid zones in the course of the simula¬ 
tion. This has a visible impact on the properties of the forward and 
reverse shocks and the associated temperature jumps, but it does 
not influence significantly the energy solved for by HERACLES. 
Very high resolution simulations that yielded better resolution of 
the shock structure gave, for example, similar results for the prop¬ 
erties of the CSM, the shock propagation speed, and the emergent 
luminosity. 

The initial configuration for the interaction is determined an¬ 
alytically. We prescribe a density and a temperature structure for 
the inner shell (tagged as the SN ejecta, but the results apply to any 
ejected shell with similar properties) and the outer shell (we assume 
the CSM arises from pre-SN mass loss in the form of a wind). We 
adopt a structure for the SN ejecta which is based on the simulations 
of core-collapse SNe we perform jPessart et al. 


, m prac- 


np: 

tice comparable to the formulation of IChevalier & Irwitil (1201 ih and 


iMoriva et alj (l2013bl) . In the present paper, we are only interested 
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Figure 2. Initial configuration for the reference interaction model X simulated with HERACLES. We show radial profiles of the velocity, the temperature, the 
mass density, and the Lagrangian mass. Global quantities and additional details for this model are given in Table[T]and Section IzT] 


in studying the basic properties of super-luminous SNe Iln and thus 
focus on one event (i.e., SN 2010jl), for which we are guided by the 
parameters of lpransson et alJ j2014h . who argue for a SN explosion 
leading to interaction with a dense and extended CSM. We do not 
include any unstable nuclei in the simulation and thus ignore any 
contribution from radioactive decay. 

For our standard interaction model, named X, we take a 
9.8 Mq inner ejecta (in homologous expansion) with 10®^ erg. Its 
density structure is given by a power law in radius with exponent 
Np =8 outside of Vo ~ 3000 km s“^ and constant within it. Its 
temperature structure is given by a power law in density with expo¬ 
nent A^t = 0.4, rising from 2000 K at the ejecta/CSM interface ra¬ 
dius located at Rt (we enforce a maximum temperature of 65000 K 
in the inner ejecta). The SN ejecta is 11.6 d old when the interaction 
starts. For the outer shell, which starts at Rt — 10^® cm, we adopt 
a wind structure with a constant velocity of 100 km s“^ (and con¬ 
stant temperature of 2000 K), but split that space into two regions 
of distinct density. Below 10^® cm, we prescribe a mass loss rate 
of 0.1 Mq yr“^, and beyond that radius we use a mass loss rate of 
10“® Mq yr~^ (i.e., two orders of magnitude smaller). The moti¬ 
vation for this is three fold. First, the CSM mass should not be un¬ 
realistically large. By truncating the high mass loss region, we can 
control the total CSM mass for any chosen mass loss. Second, the 
phase of high mass loss in massive stars like rj Car lasts for a period 
of the order of ten years, not for centuries. With our choice of outer 
radius and wind speed, the high mass loss phase has a duration of 
~ 28 yr. Third, the transition to a lower density outer region will 
facilitate a reduction in the luminosity at later times, as observed 


in SN 2010jl past 300 d (see below, and iFransson et alj|2014h . For 
convenience, we interpolate the density at the junction between the 
two shells in order to smooth the profile and avoid an overly abrupt 
change in luminosity when the shock reaches this region. The inner 
edge of the outer shell has a kinematic age of 3.2 yr. The CSM has 
a total mass of 2.89 Mq and a total kinetic energy of 5.2x 10^^ erg. 
To summarise, for the SN ejecta {R < Rt), we adopt: 

V{R) = )Rnax , 

Vo = Vmax/S , 

p{V) = poiVo/Vf'’ for V>Vo, 

p(V) = po for V < Vo and 

T{V) = 2000 (p(F)/p(Rnax))'^^ K . 

The proportionality constants po and Vmax are adjusted so that the 
total SN ejecta mass and kinetic energy match the desired values 
(see Table[T]l. For the CSM (R > Rt), we adopt: 

p(i?) = M/AtvR^V, 

V = 100km s“^ and , 

T = 2000 K 

where M = McsM.in for R < 10^® cm and M = McsM.out for 
R > 10^® cm. 

The initial conditions for model X are shown in Fig.|^ To ex¬ 
plore some dependencies of our results, we vary these ejecta/CSM 
conditions and discuss the implications in Section |5] (see also Ta- 
ble[B. 
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Table 1. Summary of simulations perfonned in this work. The reference model X is discussed in detail through the most part of the paper (Section O. while 
additional simulations (Xe3 etc.) are discussed in Section [5] For all simulations, the minimum and maximum radii of the Eulerian grid ai‘e 5x10^^ and 
2x 10^® cm, and the transition radius between the SN ejecta and the CSM lies at a radius Rt of 10^^ cm. In all cases, the mass of the inner shell (the SN 
ejecta) is 9.8 Mq, and the CSM wind velocity is 100 km s“^. Results from the HERACLES simulation for each model are given in the last three columns. The 
time to peak is the time to reach maximum from the time when the rising bolometric luminosity is only 1% of the value at peak (this way, we cancel the light 
travel time to the outer boundary, where we record the flux). We finally add results from Section[3]for the bolometric correction and colour at peak. Numbers 
in parenthesis are powers of ten. 


model 

®kin,SN 
[10®^ erg] 

Vmax.SN 
[km s“^] 

/?kin,CSM 
[10®^ erg] 

Mcsm 

[Mq] 

McSM,in McSM.out 
[M0yr-l] 

-^bol,peak 

[ergs“^] 

B.C.@peak 

[mag] 

(V^ -^)@peak 

[mag] 

^peak 

[d] 

J Ldt 
[10®^ erg] 

X 

1 

9608 

5.17(-4) 

2.89 

0.1 

0.001 

3.024(43) 

-1.06 

0.15 

19.4 

0.32 

Xe3 

3 

16642 

9.70(-4) 

2.89 

0.1 

0.001 

1.204(44) 

-1.35 

0.11 

15.7 

0.88 

Xe3m6 

3 

16642 

5.15(-3) 

17.31 

0.6 

0.006 

2.080(44) 

-1.39 

0.01 

55.7 

2.05 

Xe3m6r 

3 

16642 

6.08(-3) 

26.73 

0.6 

0.006 

1.818(44) 

-1.05 

0.06 

68.3 

2.13 

XelO 

10 

30384 

2.55(-3) 

2.89 

0.1 

0.001 

6.399(44) 

-1.46 

0.13 

12.7 

2.92 

Xel0m6 

10 

30384 

1.31(-2) 

17.31 

0.6 

0.006 

1.091(45) 

-1.80 

-0.04 

34.2 

6.89 

Xm3 

1 

9608 

1.46(-3) 

8.66 

0.3 

0.003 

3.906(43) 

-0.84 

0.13 

47.9 

0.49 

Xm6 

1 

9608 

2.87(-3) 

17.31 

0.6 

0.006 

4.751(43) 

-0.80 

0.27 

77.5 

0.63 



Figure 3. Bolometric light curve evolution for the SN Iln model X (multi- 
group; black) and its grey counterpart (blue). We add the bolometric lumi¬ 
nosity at the shock (green dashed curve), which differs at times <200 d from 
the recorded bolometric luminosity because of optical depth and light travel 
time effects (there is an 8 d light travel time to the outer boundary from the 
initial location of the interaction region and a ~ 35 d diffusion time through 
the CSM once it has become ionised). For co mparison, we show the an¬ 
alytical model light curve (red dashed curve) of iFransson et ai UlOl^ for 
SN2010jl. 


As the code is Eulerian, the flow constantly leaves the grid at 
the outer boundary. Conditions are flow out, for both the gas and 
the radiation. At the inner boundary, we use a reflecting condition 
for the radiation (zero flux). For the gas, we use an inflow condi¬ 
tion dynamically consistent with the properties of the inner shell. 
We assign to this injected material a low temperature of 2000 K, a 
velocity that preserves the homologous expansion of the inner shell 
material (V oc R), and a density pib = pix^, with x = t/{t + dt), 
and where pi is the density in the first active zone, t = Rih/Vih, 
and dt is the time step. The injected material has no dynamical or 
radiative influence on the physics of the interaction and our results. 
We continue the simulations until the interaction region moves out 
of the grid, until the luminosity becomes about 100 times smaller 
than at maximum, or when 2 years have passed, whatever comes 
first. 


2.2 Results from the Heracles simulation 

Below, we present the salient features from the interaction sim¬ 
ulation, together with illustrations of the bolometric light curve 
(Fig.®, the interaction evolution (Figs. mu, the properties of the 
photosphere and of the dense shell (Fig. [U, the evolution of the 
fluid and of the radiation properties versus radius (Fig. and the 
shock structure (fig. □oj. We define the photosphere as the location 
where the optical depth integrated inwards from the outer boundary 
is 2/3. For convenience, we consider only the electron-scattering 
opacity when computing the photosphere location. This choice is 
suitable because the bulk of the emerging radiation emerges in the 
optical spectral region where electron scattering is the dominant 
opacity source. At X-ray and UV wavelengths energy-dependent 
optical depths would need to be used. 

Initially the CSM, by construction, is at 2000 K and neutral. 
It has a huge opacity to Lyman-continuum photons but is transpar¬ 
ent at longer wavelengths (Fig. [T]l. However, after about a week, 
this CSM is ionised (the ionisation is nearly complete for hydro¬ 
gen, but only partial for helium). This rapid change occurs as an 
ionisation front sweeps the CSM (see the early migration of the 
photosphere, shown as a filled dot, in Fig. [8}, triggered by the huge 
luminosity arising at the shock. At the onset of interaction, both the 
CSM density and the ejecta velocity are maximum, so the shock 
luminosity is at its maximum. The luminosity is dominated by the 
contribution from the forward shock and can be estimated wit h 
Tshock ~ 27rr^pcsmr’sbock (see, e.g.. IChugai & Danzig 1 19941) , 
which is of the order of 3x 10'^'^ erg s“^. As time progresses, the 
fast ejecta is decelerated by the CSM and the shock luminosity 
decreases. Within 10d, the maximum velocity is ~ 6000km s“^ 
and the shock luminosity ~3x 10"^^ ergs“^. These luminosities are 
extreme. Adopting a representative radiative luminosity of a few 
10^^ erg s“^ crossing the CSM over the first ten days, this amounts 
to a total radiative energy on the order of 10"*® erg, which is 100- 
1000 times the energy required to ionise the hydrogen atoms in that 
H-rich shell. Only a small fraction of this energy is absorbed but 
this is more than sufficient to cause rapid ionisation of the CSM. 

Once the CSM is ionised, the minimum continuum opacity 
is set by electron scattering, and is moderate at all wavelengths. 
Photons, which are injected at the shock, are thermalised due to 
the huge X-ray/UV opacity and converted to UV/optical photons. 
These photons are trapped and escape after diffusing for a time 
fdiff ~ TcsmAR/c, where AR = iiphot — Reds, and ifphot is 
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Figure 4. Radial slices through the grid showing the temperature, the ratio of the internal (i.e., gas) energy and the radiation energy (i.e., the mean intensity), 
the velocity, the luminosity, the mass density, the reduced flux, the opacity (inverse mean free path) for each group (coloured cuiwes) and for electron scattering 
(dashed black curve) together with the associated optical depth (integrated inwards from the outer boundaiy; the dotted black line corresponds to r = 2/3). 
The coloured curves correspond to each one of the eight energy groups (ordered from low to high energy/frequency; see labels in bottom right panel and 
discussion in Section lzH . We use 3 groups longward of the Paschen jump, two groups in the Paschen continuum, two groups in the Balmer continuum, and 
one group for the Lyman continuum. The time is 4.1 d after the onset of the interaction. The electron scattering contribution is not included in these groups. 


the photospheric radius and Reds is the Cold-Dense-Shell (CDS) 
radius. This CDS is very narrow (i.e., 10^^ cm) and is bounded 
by the reverse shock and the forward shock. As can be seen in 
Figs.[6}{7land Fi^. | lQ| the for ward shock exhibits a Zel’dovich spike 
dZel’^vich & Raiz^ll967b where the gas temperature becomes 
very high before relaxing to its post-shock value. In practice, the 
CDS is not cold - its temperature at these early times is in excess 
of 20000 K, having, however, efficiently cooled from X-ray shock 


temperatures by radiative emission^ In our reference simulation 
(i.e., for this adopted CSM structure) at 10 days after the onset of 

^ We obtain higher spike temperatures of >10^ K with higher resolution. 
However, it is in practice difficult to resolve the spike structure given the 
scale of the problem. Shock temperatures of 10^-10^ K are obtained if we 
simulate the same interaction configuration without radiation transport (i.e., 
equivalent to hydrodynamics only). When radiation is included, radiation 
mediates the properties of the shock. The large heat capacity of the photon 
gas reduces the shock temperature. The diffusion and the escape of radiation 
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Figure 5. Same as Fig.|4] but now for a time of 29.9 d after the onset of the interaction. 


the interaction, Rphot ~ 7x 10^® cm, the CSM (electron scattering) 
optical depth is ~15, and the diffusion time from the shock to the 
photosphere is -^35 d. 

The light curve we obtain if we assume grey radiation trans¬ 
port (within the Ml model) is quite different (blue curve in Fig.[^, 
especially at early times. The luminosity is maximum at the begin¬ 
ning, because the grey opacity is very low for cold gas. In contrast, 
the opacity to X-rays and UV photons is huge at low temperature. 
Hence, in the grey case, the photons initially produced at the shock 
travel without being absorbed and reach the outer boundary unadul¬ 
terated. It is only subsequently that the streaming radiation raises 

from the shocked region exacerbates the reduction in temperature of the 
shocked region compared to a pure hydrodynamics configuration. 


the CSM temperature and increases its optical depth, causing a dip 
and a bump in the light curve, but the corresponding opacity in 
the X/UV ranges is still much lower than what it is in Nature. The 
radiation properties (bolometric luminosity, but also colours etc.) 
are thus affected. This test clearly shows that multi-group radia¬ 
tion transport (with groups that span from short to long wavelength, 
even coarsely) is important when modelling SNe Iln. 

The evolution of the photospheric radius is quite complicated. 
Initially it progresses to larger radii as the CSM becomes more 
ionised. Once the medium is fully ionised, the photospheric radius 
remains almost constant (Fig. H. Depending on circumstances, the 
photospheric radius can then either decrease or increase. First, the 
CSM may recombine and its opacity to radiation decrease, leading 
to a recession of the photosphere. This is not seen in our simulations 
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Figure 6. Same as Fig.|4] but now for a time of 179.9 d after the onset of the interaction. 


because the large and sustained supply of radiation from the shock 
maintains a high ionisation in the CSM. The other circumstance, 
which inevitably occurs, is that the shock (and the CDS), eventu¬ 
ally overtakes i?phot causing the photospheric radius to increase 
again. This occurs in this simulation at ~ 200 d, which implies that 
the average shock velocity over the hrst 200 days is ^^3000 km s“^ 
(Fig.HJ . Because of expansion, the electron-scattering optical depth 
of the CDS decreases with time.0 

The evolution of the light curve for the first ~200d is 

^ This holds once the CSM is fully ionised. Radial compression of CSM 
material by the shock does not change the CSM electron-scattering optical 
depth; this would require lateral compression, inexistent in a spherically- 
symmetric simulation. Absoiptive optical depths may increase (because 


essentially at fixed -Rphot* The morphology of the light curve 
(Fig.H differs from the monotonically decreasing shock luminos¬ 
ity because of optical depth effects, very much in t he manner of 
^^Ni powered optically- t hick ejecta of standa rd SNe l lArnettl[T982l : 
IChevalier & Irwinlboi ll : iMoriva et alJl20139l . At early times, the 
ishock is huge but the optical depth is large so radiation energy 
is stored and released on a diffusion time scale of about a month. 
Bolometric maximum occurs at f ~ tdiff. Past maximum, the SN 
radiation exceeds for a while and eventually becomes equal to the 
shock luminosity. A steady-state configuration sets in because the 
shock luminosity evolves slowly and the diffusion time is only a 

such process are often proportional to the square of the density) but except 
for the optical depth of the Lyman continuum these also decline. 
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Figure 7. Same as Fig.|4] but now for a time of 389.9 d after the onset of the interaction. 


few days (because the CDS gets closer to Tiphot)- Thus, the bolo- 
metric evolution over the first 200 d reflects primarily variations 
not in i?phot but in the radiation temperature, itself dependent on 
the amount of energy injected and stored between the photosphere 
and the CDS. For times t < fdifi, this energy accumulates and the 
temperature of the optically-thick CSM gas goes up to a maximum 
around bolometric maximum. For times t > idift, less and less 
energy is stored until the steady-steady state regime is reached. As 
the shock luminosity decreases, the radiation temperature decreases 
too and the SN luminosity ebbs0 

® In theory, the adopted value of the mean atomic weight should influence 
the gas temperature. For partial ionisation, the mean atomic weight could 
be around 0.6-0.7, rather than 1.38, which corresponds to the neutral H- 


At the end of the simulation, a total of 0.3 x 10®^ erg has been 
radiated, drawn from the 10®^ erg of kinetic energy stored in the 
inner ejecta. There still remains, untapped, 0.7 x 10®^ erg of kinetic 
energy in the system, contained in the un-shocked SN ejecta (i.e., 
the SN ejecta material not affected by the interaction, which rep¬ 
resents 50% of the total SN ejecta mass) and in the massive CDS, 
which still advances through space at ^3000km s“^. So, for this 
configuration, we obtain a conversion efficiency of kinetic to ra¬ 
diation energy of 30% (little additional kinetic energy will be ex¬ 
rich gas adopted here. We ran a simulation for the reference model X with 
a value of 0.6, but found very similar results, probably because in super- 
luminous SNe Iln, the entire material is partially ionised throughout most 
of the evolution — the ionisation state of the plasma does not change. 
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Figure 8. Illustration of the evolution of the material properties at the photosphere (left; logarithmic scale for abscissa), and at the location of maximum 
density within the dense shell (right). Note the migration of the photosphere to larger radii during the first 8 days, as well as the sudden jump in photospheric 
properties at ~ 220 d, when the CSM material exterior to the CDS is no longer optically thick (because the CDS has swept up the dense part of the CSM). 


tracted in this interaction because the remaining CSM is very low 
density). In a multi-dimensional configuration, some kinetic energy 
could also be stored in the lateral direction. Hence, the present 1-D 
simulation likely overe stimates what the equ ivalent 3-D simulation 
would produce. Unlike iMoriva et ^ j2013hh . we make no attempt 
to estimate this, in part because these lateral motions will be much 
slower than the radial motion of the material and thus cannot rep¬ 
resent a large energy loss for the radiation. 

In the next section, we will discuss in more detail the radia¬ 
tive transfer properties of this interaction. One thing to note is that 
despite the significant optical depth of the configuration, the pho¬ 
ton mean free path is non negligible. This is particularly true in the 
inner ejecta region once expansion has caused significant cooling 
because these regions do not benefit from the shock luminosity (the 
bulk of the radiation streams radially outwards from the interac¬ 
tion). There is a temperature jump ahead of the interaction (Fig. [Jo) 
caused by radiation leakage from the sh ock, as in the phenomeno n 
of shock breakout in core-collapse SNe llKlein & CheTOlieilll978h . 
This structure is comoving with the shock, since the conditions that 
cause it persist for as long as the shock remains optically-thick. 

Our results a re in agreement with the radiation hydrodynam¬ 
ics simulations of iMoriva et al.l l l2013hh. Our simulation also e m- 
phasises that the shell shocked model l Smith & McCravll2007h is 
not adequate for optically-thick sup er-luminous SNe Iln, as also 
pointed out bv iMoriva et al.l ( l2013ah . The basic inconsistency with 
the shell-shocked model is that it is not possible for a shock to cross 
a very extended CSM, ionise it, store energy within it, and subse¬ 
quently let this shocked material radiate the deposited energy. This 
scenario applies to shocks crossing the interior of a stellar enve¬ 
lope (as in successful core-collapse SNe) only because the stellar 
interior has a huge optical depth. At every location except for the 
outermost stellar layers, the radiation dominated shock progresses 
outward faster than the photons it carries in its wake because the 
photon mean free path is exceedingly small and photon diffusion 
times exceedingly long compared to the shock-crossing time. In in¬ 
teracting SNe, this configuration does not hold at all. Even for very 
high CSM densities, the shock-crossing time through the CSM will 
always far exceed the diffusion time through the shell so that radi¬ 
ation deposited by the shock is continuously being radiated by the 
shocked material. The present simulation shows that the basic light 
curve properties of super luminous SNe can be explained with a 
simple and physically consistent setup. It also illustrates how the 


simplistic “diffusion model” can be misleading, and how an argu¬ 
ment based on dimensional analysis can sometimes lead to wrong 
conclusions. 


In our simulations spectrum formation takes place primarily 
in the slow moving CSM for about 200 d, but past that the spec¬ 
trum formation region will be tied to the much faster moving CDS. 
The CDS optical depth will be too small to ensure full thermalisa- 
tion within the slow unshocked material even if the location where 
the optical depth is 2/3 (i.e., the photosphere), lies there. This im¬ 
plies that line profiles should evolve and become more Doppler- 
broadened as time passes in super-luminous SNe, while line broad- 
ening by electron s cattering would dominate only at early times 
dPessart et al.l200^ . This behaviour is important to study as it car¬ 
ries information on the potentially large residue of kinetic energy 
left untapped in the interaction. For example, this matters when es¬ 
timating the kinetic energy in the inner shell and whether this inner 
shell originates in an energetic explosion that requires the gravita¬ 
tional collapse of the iron core of a massive star. 


The radiation injected by the shock and crossing the CSM 
does work on the material and accelerates it. This radiative ac¬ 
celeration is however modest, and limited to the peak time of the 
light curve, leading to an increase of the CSM velocity from 100 
to ~200km s“^ (Fig.O. At early times, the material that receives 
a larger acceleration, just ahead of the shock, lies at a large optical 
depth and is thus less visible that the photons emitted further out, 
in regions where the acceleration is more modest. At late times, 
the material ahead of the shock is too optically-thin to receive a 
significant acceleration. The radiative acceleration in the HERA¬ 
CLES simulation is likely underestimated because of the neglect of 
line opacity and the potential boost from line desatura tion caused 
by velocity gradients jCastor et al.lll975l ; IOwoc 320151) . However, 
our simulations suggest that the radiation arising at the shock only 
cause a modest acceleration of the unshocked CSM. This is com¬ 
patible with the ne ar-constant width of the narrow P-Cygni profile 
seen in SN 2010jl jFransson et al.ll2oT^ . As we discuss in the next 
section, this modest acceleration seems an unlikely origin of the 
peak blueshift seen in the broad emission component of Ha. 
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Figure 9. Evolution of the radial profiles for the velocity, luminosity, temperature, and mass density for the reference model X. The epochs shown are 0.01,1, 
2, 3, 4, 5, 6, 7, 8, 10, 20, 30, 50, 100, 150, 200, 250, 300, 350, 400 d after the onset of the interaction. Dots refer to the location of the photosphere — such 
a location exists for about 350d, first located in the CSM and then in the CDS. Past 350d, the total electron scattering optical depth is less than 2/3, making 
electron scattering a poor frequency-redistribution mechanism at such late times. 




Radius [10^" cm] M., [Mq] 

Figure 10. Profiles for the density, temperature, velocity, and bolometric luminosity in the shock region at maximum light (all normalised for visibility) as a 
function of radius (left) and mass (right) in the reference model X. The resolution is somewhat too small to fully resolve the shock structure, a problem that 
arises from the Eulerian approach. 
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3 SPECTROSCOPIC MODELLING WITH cmfgen 
3.1 Numerical approach 

At selected epochs during the evolution, from 20 to ~ 200 d af¬ 
ter the onset of the interaction, we post-process the HERACLES 
simulation with the non-LTE radiative trans fer code CMFGEN 
jHillier & MilleJI 19981: iHUlier & Dessarl2012h . While there is no 
explicit time-dependence in the CMFGEN simulations undertaken 
here, time dependence is implicitly taken into account since we 
use the density, temperature, and velocity structures computed with 
H ERACLES. We t hus re tain a much higher physical consistency than 
in IPessart et al.l ( l2009h . We use an adaptive grid to resolve care¬ 
fully the strongly varying density and temperature profiles, as well 
as the optical depth. The grid is truncated at small radii when the 
inward integrated Rosseland-mean optical depth reaches 30. The fi¬ 
nal models discussed here use 100 grid points and adopt a turbulent 
velocity of 20 km s“ ^. 

We use the same H, He, and Fe mass fractions as for the 
radiation-hydrodynamical simulations. We limit the model atom to 
H I, He I, He H, Fe I-V I with the same super-level assignment as 
in IPessart et al.l ( l2013h . Higher ionisation stages are unnecessary 
since the gas temperature in the simulations at the times considered 
is below ~ 80000 K (this maximum temperature is limited to a few 
zones at the shock and to early times). In the present simulations, 
the spectrum formation region has a characteristic temperature of 
about 10000 K or less. A discussion of the signatures of ions from 
CNO and other intermediate mass elements (which are merely trace 
elements compared to H and He here) is deferred to a subsequent 
study. 

The current version of CMEGEN does not allow us to calcu¬ 
late non-LTE line blanketed model atmospheres for non-monotonic 
velocity flows. This property of interacting SNe is however funda¬ 
mental and should be taken into account. In this work, to overcome 
this limitation and take into account the non-monotonic velocity 
structure, we ignore lines when treating the continuous spectrum, 
but model the lines for the non-LTE analysis using the Sobolev 
approximation. Non-monotonic velocity flows do not cause major 
issues with the Sobolev approximatiorlj since it is a local approxi¬ 
mation that uses the a bsolute value of the velocity gradient to com¬ 
pute the optical depth ( ICastoill970l) . Since the density and temper¬ 
ature structures are taken from the hydrodynamical simulation we 
believe this to be a very reasonable first approach. 

To compute the observed spectrum we wrote a new routine ca¬ 
pable of computing J in the comoving frame in the presence of non¬ 
monotonic velocity fields. With this routine we were able to allow 
for the influence of (non-coherent) electron scattering which was 
iterated to convergence using a simple lambda-it eration. Relatively 
minor changes were then made to CMF_FLUX teusche & Hillierl 


® One problem we needed to address is that the Soholev optical depth be¬ 
comes infinite for one particular direction when the radial velocity gradient 
is not positive. This occurs because in a spherical flow the Sobolev optical 
depth is proportional to 1/|1 + cr/r^|, where a = d\nV/ din R — 1. The 
angular factor will be infinite when fi = ibl/\/—cr. In practical terms 
this does not affect the escape probability very much due to the small 
solid angle subtended by the singularity. However, in the standard com¬ 
putational approach the angle grid is not well sampled (the sampling be¬ 
comes worse as we move away from fj, = 1), and numerical instabilities 
occurred at one or two depths, which we overcame by limiting the increase 
in optical depth. In this work we used the simple approach of choosing 
max(|l + cr/r^l, —0.2(1 + a)) for |1 + cr/r^| when cr < 0. 


l2005h to facilitate the computation of the observed spectrum in the 
observe r’s frame. _ 

In IPessart et alJ (l2009li . we used a prescription for the den¬ 
sity structure of the emitting region and the ansatz was that an 
optically-thick shell was radiating the flux. In that approach, the 
entire grid was in radiative equilibrium, and the flux was forced 
to diffuse through the densest regions of the atmosphere (these re¬ 
gions were also the hottest by imposing the diffusion approxima¬ 
tion at the base). The dynamical configuration presented in Sec- 
tionl^is obviously quite different, but this simplistic treatment is in 
fact not so bad at early times, when the interaction region and the 
inner CSM are optically thick because in this case the details of the 
inner boundary conditions are lost. Hence, f or as long as the con- 
figuration is optically thick, the approach of IPessart et al.l (l2009h 
can provide a good match to the observations of some SNe Iln like 
1994W without being dynamically consistent. 

3.2 Spectral evolution 

We thus compute non-LTE spectra at 20.8, 41.7, 62.5, 83.3, 104.2, 
125.0, 145.8, 166.7, 208.3 d after the onset of the interaction. The 
spectral evolution shown in Eig. reveals a striking similar¬ 
ity to observat ions of numerous super-luminous SNe Iln, includ¬ 
ing SN 2010jl dZhang et al.ll2012l : lFransson et alj|2014l) . The SEP 
evolves very slowly with time, with a progression towards redder 
colours, a diminishing flux shortward of 4000 A, and a decreasing 
ionisation visible through the strengthening of Ee II lines. H I and 
He I lines behave in a more complicated way. 

Line profiles prior to 100 days show the distinct symmetric 
morphology resulting f rom frequency redistribution by sc attering 
with thermal electrons dChugailboOll : IPessart et al.ll2009l) . with a 
typical width of a few 1000 km s“^ — in earlier calculations, where 
we limited the velocities to 100km s“^ on the CMFGEN grid, we 
obtained similar profiles. The width of this component only reflects 
redistribution of photons in frequency space by electron scattering 
- it provides no information about the shock velocity. Superposed 
on the broad wings, is a very narrow component with a P-Cygni 
profile morphology. This narrow component corresponds to line 
emission and absorption weakly affected by electron scattering, 
and which preferentially arises from the outer regions of the CSM 
(which are moving at ^100km s“^ in our model X) where line 
opacity may still be large but where the electron-scattering opacity 
is negligible. 

At a qualitative level, the early-tim e spectra shown in Fi g.llll 
are analogous to those computed in IPessart et al.l d2009l) de¬ 
spite the very different approaches in the description of the den¬ 
sity/temperature structure. 

After 100 days, the line profile morphology changes. In addi¬ 
tion to the two components discussed above, we see a blue shifted 
emission component that arises from the CPS. It is blue shifted be¬ 
cause the CPS is moving at a velocity of ~3000km s“^ and is op¬ 
tically thick in the continuum - line photons emitted in the red are 
absorbed by the near side of the CPS. T he blueshift of peak emis- 
sion seen in non-interacting Type II SNe jPessart & Hillieil2005bl) 
has a similar origin. These Ha properties agre e well at the quali - 
tative level with the observations of SN 2010jl dZhang et alj|2012h . 
although the magnitude of the emission blueshift is overestimated 
and the magnitude of the line width is underestimated (Fig. □3. 

We give a summary of the CMFGEN model properties in Ta- 
ble|2 The CMEGEN bolometric luminosity agrees with the HERA¬ 
CLES value to within 5-25% at early times (i.e. prior to ~200d), 
which is noteworthy since the two codes differ vastly in the treat- 
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Figure 11. Montage of spectra at selected epochs (see right label) for our reference model X. The high-frequency patterns in the UV are caused by line 
blanketing from Fell and Felll in this model, while the optical is composed of Hi, He I, and Fell lines. At early times, lines typically show a hybrid 
morphology with a narrow core and extended wings, up to about 2000 km s~^ from line centre. These wings are due to incoherent electron scattering. As time 
progresses, a blue shifted component, associated with emission from the CDS, grows in strength. Line broadening is then caused by electron scattering and 
expansion. 
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Table 2. Summary of properties for the CMFGEN simulations based on the reference interaction model X. We include the properties at the photosphere, at the 
location of maximum density in the interaction region (the CDS), and the properties of the fitted blackbody (see Fig. ll3l for an example). We select epochs 
when the CDS is well below the photosphere, so that the interaction region is optically thick. 


age 

-^bol 

Mb My Mb Mj 

-^phot 

T'phot 

-^e,phot 

-^cds 

Teds 

-^e,cds 

^bb 

Thb 

Thb 

[d] 

[ergs-l] 

[mag] [mag] [mag] [mag] 

[10l®cm] 

[10^ K] 

[cm“®] 

[10^® cm] 

[10^ K] 

[cm-3] 

[10l®cm] 

[lO'i K] 

[ergs-l] 

20.8 

2.15(43) 

-18.61 -18.68 -18.85 -18.86 

7.17 

1.33 

4.36(8) 

1.25 

1.66 

1.76(11) 

1.172 

1.084 

2.259(43) 

41.7 

2.23(43) 

-18.84 -18.98 -19.14 -19.21 

7.10 

1.27 

4.42(8) 

1.78 

1.18 

1.33(11) 

1.448 

1.033 

2.371(43) 

62.5 

1.21(43) 

-18.40 -18.68 -18.90 -19.03 

7.02 

1.09 

4.05(8) 

2.86 

1.01 

1.05(11) 

1.616 

0.879 

1.289(43) 

83.3 

9.88(42) 

-18.16 -18.55 -18.82 -19.01 

7.04 

1.00 

4.00(8) 

3.37 

0.79 

3.93(10) 

1.848 

0.788 

1.050(43) 

104.2 

8.65(42) 

-17.95 -18.43 -18.74 -18.99 

7.11 

0.97 

3.94(8) 

3.91 

0.76 

2.25(10) 

2.083 

0.715 

9.068(42) 

125.0 

7.99(42) 

-17.74 -18.27 -18.62 -18.89 

7.08 

0.98 

3.96(8) 

3.64 

0.61 

9.24(9) 

2.199 

0.673 

8.496(42) 

145.8 

7.26(42) 

-17.56 -18.12 -18.49 -18.77 

7.15 

1.00 

3.90(8) 

4.94 

0.69 

9.51(9) 

2.199 

0.646 

7.503(42) 




Figure 12. Montage of multi-epoch spectra for Model X (left) and SN2010jl fright: l^ang et al1l2012h shown in velocity space with respect to the rest 
wavelength of Ha. Qualitatively, the model reproduces well the strengthening of the blue-shifted emission component and the increasing width of the observed 
Ha line profile. Quantitatively, the magnitude of the blueshift is somewhat too high while the width of the emission feature is underestimated. 
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Figure 13. Comparison of the total synthetic flux (black), the contin¬ 
uum synthetic flux (blue) and the best match blackbody flux to the optical 
range (red), for the reference model X at 41.7 d after the onset of inter¬ 
action. In the same order, the cumulative luminosities are 5.83x 10® Lq, 
6.18x 10® Lq, and 4.45X 10® Lq. The blackbody radius and temperature 
are 1.45x 10^® cm and 10331.6K. The photospheric radius and tempera¬ 
ture at that time are 7. lOx lO^^cmand 12670.0 K. 

merit of the radiative transfer (Ml method using eight energy 
groups versus a two-moment solver using 10® energy groups) and 
of the thermodynamic state of the gas (LTE versus non-LTE) — 
past 200 d, the holometric luminosity predicted by the two codes 
differs sizeably (about a factor of 2), probably because non-LTE 
effects are stronger (the CDS is no longer embedded within the 
optically thick CSM). A mismatch is expected during the first 1-2 
months because light-travel time effects are ignored in the compu¬ 
tation of the flux in CMFGEN (the time dependence of the model 
is implicitly taken into account since we fix the temperature, but 
no explicit time dependence is accounted for when computing the 
flux in the formal solution of the radiative transfer equation in CM¬ 
FGEN). This matters when the properties of the interaction vary on 
a short time scale (i.e., comparable to the free-flight time of ~8d 
to the outer boundary, where the luminosity is recorded). 

The slow evolution of the SED reflects the slowly varying 
shock luminosity (which arises from the smooth variation in both 
the ejecta density/velocity profile and from the slowly varying CSM 
density for a wind configuration), leading to a decrease in radi¬ 
ation/gas temperature in the optically thick regions of the CSM. 
Consequently, the electron scattering wings weaken as the ionisa¬ 
tion level and the electron scattering optical depth decrease (this is 
caused by the shrinking of the spectrum formation region, which is 
bounded by the fixed photosphere and the outward moving dense 
shell or shock). For about 200 d after the onset of the interaction, the 
spectrum formation process is essentially unchanged in this simu¬ 
lation. 

It is customary in the community to compare the SED with a 
blackbody. At each epoch, we fit the CMEGEN SED flux ratio be¬ 
tween 4800 and 8000 A to infer Thh — we use an iterative proce¬ 
dure to find the temperature for which the model and the blackbody 
have between these two wavelengths the same flux ratio to within 
one part in 10000). We determine the corresponding blackbody ra¬ 
dius iZbb by matching the blackbody to the model luminosity L\ at 
A (i.e., through i?bb = L\/B\{T’o’a)). The total blackbody lu¬ 
minosity Lbb is then An example is shown in Fig. 1131 

Although the match is excellent in the fitted range, the SED is not a 



Figure 14. Radial variation of the optical depth associated with electron 
scattering, various continua, and lines of Ho and He 15875 A (for which we 
plot the Sobolev optical depth). Thermalisation in the Lyman continuum is 
ensured, but not in the Balmer and Paschen continua. The large Sobolev 
length caused by the small velocity gradient in the CSM causes the Ho line 
optical depth to be very large, and to exceed that in the Lyman continuum. 


blackbody. We find non-planckian effects here associated with di- 
lution due to electron sc attering (causing Rhh < i^phot; see, e.g., 
IPessart & Hillieill2005ah . and flux excess rather than blanketing in 
the UV (as if multiple blackbody emitters with different tempera¬ 
tures were contributing). 

We show the optical depth for electron scattering, selected 
continuum wavelengths, as well as selected lines in Fig. [14] It is ev¬ 
ident that thermalisation is complete for UV photons in this model, 
but photons redder than the Lyman edge poorly thermalise, even at 
an early epoch of 41.7 d after the onset of interaction. For these, 
electron scattering is the main source of opacity. The blackbody 
radius underestimates the radius of the CDS at all times, although 
the bulk of the SN flux arises from outside the CDS. Overall, the 
radiation field in this SN configuration is not well described by the 
Planck distribution. These considerations are confirmed by Fig. [B] 
which illustrates how much smaller the mean intensity is compared 
to the Planck function, except at the location of lines. Furthermore, 
in our reference model, the bulk of the flux arises at all times from 
layers where the electron scattering optical depth is never greater 
than 10, and more typically a few only (see also Fig. do). 

The shallow density structure, the wavelength-dependent 
opacities, and the strong albedo, lead to a very spatially extended 
spectrum formation region. For example, the emission/absorption 
of Ha photons occurs throughout the grid, with the exception of 
the outermost regions of the CSM where the density is very low. 
At early times observed line photons primarily originate at larger 
radii because line photons emitted deeper suffer more scatterings, 
and hence have a greater probability of being absorbed. Line pho¬ 
tons that could accrue very large frequency shifts by multiple scat¬ 
tering with thermal electrons tend to be absorbed before they can 
escape. In our reference model, we find that the bulk of Ha photons 
arise from regions with an electron scattering optical depth of a few, 
giving rise to emission in line wings extending to ~ 3000 km s“^ 
from line centre, and characterised by a profile with a FWHM of 
^1500km s~^. T his is somewhat similar to values in SN2010jl 
dZhang et alj2012l) at early epochs. However at later times the pro¬ 
files in SN 2010jl are broader and blue shifted. We also see this 
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Figure 15. Left: Variation of the ratio of the mean intensity Jx and the Planck function Bx in the simulation at 41.7 d, shown as a function of wavelength and 
ejecta/CSM location. The label at right gives the electron-scattering optical depth res associated with each panel. At this time, the CDS is at res ~ 10. Right: 
Variation of Jx and Bx versus radius and shown for various continuum wavelengths. 


effect at late times, with the extra emission associated with the 
CDS (Fig. ll2t . The blueshift of the broa d Ha line is caused by an 
optical depth effect t hat affects all SNe jPessart & Hillieill2005bl : 
lAnderson et al.l2Ql^ . At early times, the peak blueshift is not seen 
because the Ha photons emitted from the CDS are unable to es¬ 
cape - they undergo numerous electron scatterings and are eventu¬ 
ally destroyed by a continuum absorption. All these properties are 
illustrated in Figs. fToHTT] 

As found in earlier studies (e.g.. lHillieill99lh the adopted tur¬ 
bulent velocity also affects the strength of the electron scattering 
wings. This occurs since a larger turbulent velocity increases the 
Sobolev length, and hence leads to an increased probability that a 
line photon scattering within the Sobolev resonance zone will ex¬ 
perience an electron scattering interaction, which, because of the 
large Doppler shift, will cause it to be shifted out of resonance al¬ 
lowing it to escape. 


4 POLARISATION CALCULATIONS 

The dense shell that forms in the interaction is Rayleigh-Taylor 
unstable and should be simulated in 3-D — this is left for future 
work. There is, however, observational evidence from spectropo- 
larimetric observations that some SNe Iln are asymmetric on large 


scale I Leonard et alj|20o0: IWang et ^l200ll : iHoffman et al.|[200^ : 


IPatat et alj201lh . This we try to address in this section by perform¬ 

ing 2-D polarisation calculations. 


We post-process the CMFGEN simulation for model X at 41.7 d 
by compu ting linear polarisation pr ofiles in the same fashion as de¬ 
scribed in lPessart & Hilli^ ( 1201 ih . We break spherical symmetry 
by enforcing a latitudinal variation in density (thus preserving axial 
symmetry). In practice, we use p(r, /r) = p(r, /r = 0)(1 + Aipd), 
where p = cos 9, and 6 is the colatitude. The magnitude of the 
asymmetry is controlled through the parameter Ai, which we vary 
from 0.1, to 0.2, 0.4, 0.8, and 1.6. The resulting axially symmet¬ 
ric ejecta has therefore a prolate morphology, with pole-to-equator 
density ratios (at a given radius) of 1.1, 1.2, 1.4, 1.8, and 2.6. This 
adopted type of asymmetry is chosen for its simplicity - other 
choices might also yield similar continuum and line polarisations. 


We show the flux in the Ha line and in the overlapping 
continuum in Fig. [T^ for an observer along a line-of-sight per¬ 
pendicular to the axis of symmetry. Along this direction, the 
continuum polarisation has its maximum value, reaching about 
1.6% for Ai = 1.6. As we progress towards the Ha rest wave¬ 
length from continuum regions (i.e., > 3000 km from line cen- 
tre), the polarisat i on decreases steadily , as observed in SN2010jl 
dPatat et al.ll201 ll : IWilliams et al.ll201^ and in some Wolf-Rayet 
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Figure 16. left: Illustration of the wavelength (A) and depth {R) dependence of the quantity A), where 5F{R,X) = 

(27r/Z)^) As: 77(p, 2 :, A) and i?i 5 is R in units of 10^^ cm. Here, D is the distance; A/? and A 2 : are the shell thickness in the ra¬ 

dial direction and along the ray with impact parameter p, respectively; p(p, z, A) and r(p, 2 , A) are the emissivity and the ray optical depth at location (p, z) 
and wavelength A. The gray scale in the bottom panel shows how varies as we progress inwards from the outer boundary, indicating 

the relative flux contributions of different regions. If we choose i?i5 as the minimum radius on the CMFGEN grid, we recover the total flux. The grayscale 
also shows how lines increasingly broaden (and the peak line flux decreases) as we progress outwards. The dash-dotted line corresponds to the radius of the 
electron-scattering photosphere, and the dashed line (which overlaps with x-abscissa here) corresponds to the CDS radius. The top panel shows selected cuts 
(see right label) of the quantity ^)- From top to bottom (/2i5 = 10, 7.2, 5.0, 3.0), A) represents 4, 45, 77, and 98 per cent 

of the total emergent flux at 6800 A. Right: Same as left, but now zooming in on the Ha region. The depth variation of the line width, influenced by electron 
scattering, is evident in the bottom panel. 


stars dSchulte-Ladbeck et al .ll99lHHarries etal.ll998h . In our sim¬ 
ulation, this reduced polarisation across the line is a signature of 
the distinct formation of the continuum and the line. Because of the 
lower optical depth, continuum photon arise from a much greater 
electron-scattering optical depth than do Ha line photons (Fig. 114b . 
Furthermore, line photons emitted in a recombination process start 
off with zero polarisation, and therefore can only accumulate a 
lower polarisation through a reduced number of scatterings before 
escape. This depolarisation is maximum within 100 km s“^ of line 
centre since over this spectral range, one sees photons that under¬ 
went essentially no scattering with free electrons (photons scattered 
by free electrons receive a wavelength-shift and, statistically, en d 
up somewhere in the wings; see discussion in lPessart et ahlEoP^ . 
In this context, whatever the imposed deformation on the ejecta, the 
polarisation at line centre is zeroQ 

Surprisingly the polarisation is parallel to the major axis (i.e., 
positive; see top panel of Fig. d — this property is, naively, not 


^ Equivalent simulations for Type II SNe may yield a residual polarisa¬ 
tion at line centre because thes e photons are emitted a t a non negligible 
electron-scattering optical depth IPessart & HillieJlOl ih . Spectrum forma¬ 
tion in SNe II-P and SNe Iln is drastically different. 


expected since scattering of the light from a point source by an opti¬ 
cally thin pr olate spheroid has the po larisation at right angles to the 
m ajor axis jBrown & McLeanl ll977h . In our convention (and that 
of iBrown & McLe^ 1977 ), the polarisation is positive when it is 
parallel to the axis of sy mmetry (for further details, see Section 2 of 
IPessart & HilliedlTOllll . The reason for this behaviour in the cur¬ 
rent simulation is that the CSM is optically thick, and the observed 
(and hence scattered) flux predominantly escapes from equatorial 
regions, where the densities are lowest (Fig|19t. 


The correspondence between our simulations and the mor¬ 
phology of the polarised and of the total fluxes of SN 2010jl sug¬ 
gests first that our model for super-luminous SNe Iln like 2010jl 
is adequate. Secondly, the observed polarisation implies that the 
ionised CSM causing the polarisation is significantly asymmetric. 
If axially symmetric, a pole to equator density ratio of ~ 3 could 
explain the observed continuum polarisation. While the CSM is 
asymmetric, it is however unclear whether the inner shell is also 
asymmetric — it does not need to be. Multi-dimensional radiation 
hydrodynamics simulations are needed to model the spatial prop¬ 
erties of the shocked CSM (for example, in the current context, the 
shock strength should also be latitude dependent, causing variations 
in temperature, ionisation, velocities etc with angle), so that we can 
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Figure 17. Same as Fig. 1161 but now for a time of 208.3 d after the onset of interaction. In the left column, the greyscale is shown on a logarithmic scale. From 
top to bottom (fJi 5 = 10, 7.2, 6.4), SF(R, A) represents 6, 41, and 94 per cent of the total emergent flux at 6800 A. The flux contributions from 

the fast moving CDS and the slow moving CSM are clearly visible. The emission from the CDS is blue shifted, while the CSM emission is symmetric around 
line centre. Electron scattering affects the emitted photons from both components. 


better characterise the nature of the asymmetry at the origin of the 
observed polarisation in SNe Iln like 2010jl. 


5 DEPENDENCY ON SOME MODEL PARAMETERS 

Varying the properties of either the CSM mass or the explosion 
energy can affect the emergent radiation. This has been discussed 
in t he past for the bolometric luminosity or the colours (see , 

e.g..lChugai et al.l2004l^oslev et alj200 ^•. Ivan Marie et alj2010l: 

iMoriva et alJ^OlS^T Here, we repeat such explorations to investi¬ 

gate the impact on spectral signatures as well. 

5.1 Impact on the bolometric light curve 

Taking model X as a reference, we produce models in which the 
SN ejecta kinetic energy is scaled by a factor of 3 (Xe3) or 10 
(XelO). We produce models in which the CSM density is scaled by 
a factor of 3 (Xm3) or 6 (Xm6). We also consider combinations of 
both, with models Xe3m6 or XelOmb. Finally, we consider model 
Xe3m6r, where the high-density CSM extends to l.SxlO’^^cm 
rather than 10^® cm only. Table [T] summarises these model prop¬ 
erties. 

These alterations produce a considerable diversity of bolomet¬ 
ric light curves (Fig.1201. in directions that are well understood. In¬ 
creasing the kinetic energy causes a much greater shock luminosity 
at all times, but since the CSM structure is not changed, the diffu¬ 
sion time through the CSM is unchanged (the time to maximum is 
decreased by just a few days for higher SN ejecta kinetic energy). 


The bolometric light curves shift upward for increasing energy. 
The mean radiation-energy density in the CSM also increases with 
shock luminosity, which raises the radiation temperature. Colours 
are thus strongly affected, getting bluer for higher energy (see be¬ 
low and Table [T] where we give the bolometric correction and the 
colour V — 1 al bolometric maximum). 

Varying the CSM mass (corresponding in our approach to a 
change in mass loss rate) modulates the impedance of the “interac¬ 
tion engine”. For larger CSM density and total mass, more kinetic 
energy is extracted from the SN ejecta. This energy is converted 
primarily into radiation. Consequently, varying the CSM mass af¬ 
fects the light curve properties. Higher CSM mass implies higher 
CSM optical depth and a longer diffusion time, hence a delayed 
time of peak. A larger CSM mass buffer implies more kinetic en¬ 
ergy is tapped so the bolometric maximum is more luminous. In the 
model set X, Xm3, Xm6, the total energy radiated is 0.32, 0.49, and 
0.63 X10®^ erg (for the same ejecta kinetic energy of 10®^ erg), so, 
the conversion efficiency can be considerably enhanced when the 
CSM mass exceeds significantly the SN ejecta mass, with values 
here as high as ~ 70% (models Xe3m6, XelOmb, or Xm6). 

Combinations of higher CSM mass and higher ejecta ki¬ 
netic energy (e.g., model Xe3m6) ) shift the reference model 
(quite suitable to S N 2010jl) to the domain where SN 2006gy lies 
I Smith et aHl2007h . with peak luminosities of the order of a few 
10^^ erg s“^ sustained for 2-3 months and a peak time that occurs 
later (~2 months after the onset of the interaction). 

Finally, extending the region of high density in the CSM from 
10^® to 1.5x10^® cm (models Xe3m6r versus model Xe3m6) leads 
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Figure 19. Colour map of the specific intensity (scaled by the impact pa¬ 
rameter p) on the plane of the sky at —3992 km s~ ^ from the Ho; line centre 
(this Doppler velocity corresponds to a continuum region — see top panel 
of Fig.|T§), for the model with Ay = 1.6, at 41.7d after the onset of the 
interaction, and seen at a 90° inclination (i.e. edge-on). We overplot seg¬ 
ments (black bars) showing the local polarisation strength and orientation 
(the length corresponding to 100 per cent polarisation is shown at the bot¬ 
tom right). 


Figure 18. Top: Total and polarised flux across the Ha line and the over¬ 
lapping continuum for different inclinations (90° corresponds to an equa¬ 
torial view), at 41.7 d after the onset of the interaction (hence close to the 
epoch of maximum light), showing the larger continuum polarisation, the 
depolarisation across the line profile, and the unpolarised line core. The 
prolate asymmetry is characterised by a pole-to-equator density ratio of 2.6 
(Ai = 1.6). Middle: Same as top, but now showing the percentage polar¬ 
isation for different magnitudes of the asymmetry {A\ values between 0.1 
and 0.6) and for an equatorial view. Bottom: Spectro-pol arimetric observa - 
tions of the Ha line for SN 2010jl on November 18.2UT jPatat et al.l201lh . 
which is within a few days of our model time (left panel). In the context of 
our model, the observed polarisation may be interpreted as arising from a 
prolate CSM with a pole-to-equator density ratio of ~ 3. 


to a sustained high luminosity beyond 300 d. In fact, all the light 
curves in Fig 1201 show a break after 150-300 d, except for model 
Xe3m6r. This light curve break in our simulations is associated 
with the shock (or the CDS) reaching the outer edge of the denser 
parts of the CSM (the shock in model Xe3m6r has not yet reached 
1.5x 10^® cm at the last time shown in Fig. 120b and starting to en¬ 
counter lower density material (see Table^for details). 

The pre-acceleration of the CSM by the radiation from the 
shock depends on the adopted configuration. In model X, we 
obtained a modest enhancement in CSM velocity from 100 to 
200 km s“^. This pre-acceleration is stronger when we use a larger 
ejecta kinetic energy (stronger shock, higher luminosity) or if the 
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CSM density is increased (there is a delicate balance to get here 
since a higher CSM density means a higher opacity and thus more 
radiation momentum extracted but it also implies more inertia) — 
see Fie.l2n 


5.2 Impact on maximum light spectra 

The diversity in bolometric light curves discussed in the preced¬ 
ing section is also strongly present in the spectra at maximum light 
(Fie. \22l . For the highest energy explosion, the shift in temperature 
in the ionised CSM leads to a major rise in ionisation for hydro¬ 
gen and helium, producing very blue and nearly featureless spectra 
(whatever blanketing from Fe occurs in the UV and far-UV). For 
models with enhanced CSM mass, the temperature/ionisation con¬ 
ditions are comparable to reference model X, but the line fluxes are 
increased. The higher the CSM mass, the redder the colour at peak. 

Interestingly, in models with very high CSM mass, although 
the B aimer line fluxes are much higher, the extent of the elec¬ 
tron scattering wings does not increase beyond a maximum of 
~ 3000 km s“^. The likely reason is that line photons do not just 
scatter off free electrons in the CSM — there is a non-zero prob¬ 
ability that they are absorbed in a bound-free or free-free process 
leading to photon destruction. The broader and blue shifted lines 
seen at late times require that we see photons blue shifted by ex¬ 
pansion. For the present case, the broader blue-shifted component 
seen at late times arises from the CDS. We see a blue-shift since the 
CDS absorbs emission from the far side of the shell. Thus the Ha 
blueshift seen in SN2010jl iSmith et al.ll2012l : IZhang et al.ll2012l : 
iFransson et alj20l4h might simply just stem from an optical-depth 
effect, without any dust intervening in the process. As the optical 
depth of the shell declines, we expect the profile to become more 
symmetric. 


6 CONCLUSION 

In this work, we have presented numerical simulations for inter¬ 
acting SNe in which both the kinetic energy of the explosively- 
produced inner shell and the mass of the wind-produced outer shell 
are very large. We used an inner shell kinetic energy of 10®^ erg or 
more, and a CSM mass of at least 3 M©. These simulations corre¬ 
spond to possible configurations for super-luminous SNe Iln. 

Our contribution is novel because we post-process our multi¬ 
group radiation-hydrodynamics simulations with radiative trans¬ 
fer tools, to compute non-LTE spectra and polarisation signatures. 
With this more global approach, we can study the numerous fea¬ 
tures that make super-luminous SNe Iln unique. 

We find that the light curve of a super-luminous SN Iln is com¬ 
posed of several distinct phases. Initially, the cold CSM is neutral 
and transparent to optical photons. As the interaction starts, the ra¬ 
diation from the shock drives an ionisation front through the CSM 
and makes it optically thick at all wavelengths. This first phase 
takes about a week after the onset of the interaction, which cor¬ 
responds roughly to the light crossing time through the shell. At 
the time we record the first emerging photons from the SN Iln, the 
CSM is therefore already optically thick. In our reference model, 
the photosphere is then located at 7x 10^® cm and will remain at 
that radius for months. The second phase corresponds to the shock 
crossing of the optically thick CSM. In our reference model, this 
phases lasts > 200 d. Over this duration, the shock is located below 
the photosphere, at an electron-scattering optical depth of about 
15 initially. Radiation from the shock is thus continuously released 


deep in the CSM and comes out on a diffusion time scale. This con¬ 
text is analogous to optically thick ejecta influenced by radioactive 
decay, although the shock energy is not all channeled into radiation 
energy. Here, the shock luminosity decreases with time because 
of the decreasing shock strength as more and more CSM is being 
swept up. Although the CSM hardly expands, the diffusion time 
scale and the optical depth decrease because of the expansion of the 
interaction region. We thus obtain a bell shape morphology for the 
light curve around maximum. This variation is not associated by the 
migration of the photospheric radius, but by the variation in shock 
luminosity, radiation temperature, and CSM/photospheric temper¬ 
ature. Optical depth effects on the light curve persist for as long as 
the shock is located below the ~ 7x 10^® cm. The third phase starts 
when the shock overtakes this location. The photosphere then fol¬ 
lows the dense shell. 


Our multi-group radiation hydrodynamics simulation pro¬ 
duces the basic light curve morphology of type Iln super-luminous 
SNe. Our reference mode l was tuned to roughly match the inferred 
properties of SN2010jl l lFransson et alj|2014h . but other simula¬ 
tions with different inner-shell kinetic energy and CSM mass en¬ 
compass a variety of light cu rves, including extreme events like 
SN2006gy l lSmith et alj2007l) . 


We illustrate the importance of performing multi-group radi¬ 
ation hydrodynamics. By contrast, grey (i.e., one-group) radiation 
hydrodynamics simulations do not account adequately for the opac¬ 
ity of the material and produce a discrepant SN Iln light curve mor¬ 
phology. 


Ou r simulations confirm the earli er work of iMoriva et'H] 
( l2013bh that the shell shocked model of I Smith & McCray l l2007h 
is not physically sound. The reason is that in interacting SNe, eyen 
a massive CSM cannot have a continuum optical depth greater than 
a few tens. Because such CSM are by nature extended, the shock¬ 
crossing time far exceeds the time it takes the shock luminosity 
to diffuse through the optically thick layers. The situation is anal¬ 
ogous to shock emergence in an exploding star, except that here 
the configuration is quasi steady state — it persists for as long 
as there is optically-thick CSM material ahead of the shock. Our 
re sults, which a r e in ag reement with the independent simulations 
of iMoriva et alj l l2013bh . provide a simple and physically consis¬ 
tent explanation for the light curves of type Iln super-luminous 
SNe. Overall, the shell-shocked model really belongs to the domain 
of (successful) core-collapse SNe in which a radiation-dominated 
shock crosses a stellar envelope whose optical depth is huge every¬ 
where apart from the surface. 


The non-LTE radiative transfer simulations, based on the ra¬ 
diation hydrodynamics simulations, produce spectra that are ini¬ 
tially blue and redden progressively as the CSM/photosphere cools 
down after bolometric maximum. We obtain symmetric emission 
profiles, broadened by non-coherent scattering with free electrons 
(of thermal origin), with wings that extend up to 3000 km s“^ 
from line centre, and together with a narrower P-Cygni profile 
(< 100km s“^) forming in the outer cold CSM. Varying the ejecta 
kinetic energy or the CSM mass alters the shock luminosity, the 
mean radiation energy in the optically thick CSM, and causes clear 
variations in colour and ionisation. Varying the CSM density by a 
factor of 6 changes the strength of recombination lines but does 
not produce broader lines at bolometric maximum, despite the en¬ 
hanced electron-scattering optical depth. Photons emitted at an op¬ 
tical depth r perform scatterings before escaping. If the albedo 
is 99%, photons coming from r > 10 are likely absorbed within the 
ionised CSM and photons do not obtain arbitrarily large frequency 
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Figure 20. Light curve diversity obtained by varying some of the parameters of our reference model X (see details in Table[T). Models whose name contains 
en have a SN ejecta kinetic energy scaled by a factor n relative to model X. Models whose name contains mn have a CSM mass loss rate scaled by a factor n 
relative to model X. Model Xe3m6r differs from the rest because the dense part of the CSM extends to 1.5x 10^® cm rather than 10^® cm. 


shifts from line centre. Broader electron scattering wings may arise 
for a greater albedo. 

Our simulations indicate that photon absorption is probably 
efficient only in the Lyman continuum. At low photon frequen¬ 
cies, absorptive processes are weak and therefore photon absorp¬ 
tion occurs only if the total optical depth at the site of emission is 
large (say above 10) — in our reference simulation with a mass 
loss rate of 0.1MQyr“^ for the CSM, the maximum CSM con¬ 
tinuum T is only 15. Early on, the optical depth in the CDS is high 
enough to ensure complete thermalisation. Then, line broadening is 
probably dominated by non-coherent scattering with thermal elec¬ 
trons. However, at later times, thermalisation does not occur at all 
depths in the CSM. In the reference model, this is associated with 
an increased broadening and blueward asymmetry of line profiles 
with time. These features are also present in the observations of 
SN2010jl (e.g., Ha; Fig. [Hi — this suggests a growing contri¬ 
bution from Doppler broadening, connected with photon emission 
from the fast expanding CDS. We believe this hybrid line forma¬ 
tion is at the origin of the very large Ha line width in SN 2010jl, 
as well as the blueshift of its peak emission. Interestingly, the mod¬ 
est peak luminosities and the absence of broad lines at all times 


in events like SN 1994W jChugai et al.| l2004 IPessart et al.ll200^ 
or SN 201 Iht jHumphrevs et al.l 2012h are perplexing because the 
emission from the CDS should be seen once the CSM becomes 
transparent. 

To complement our radiation hydrodynamics simulations with 
HERACLES and our non-LTE radiative transfer simulations with 
CMEGEN, which assume spherical symmetry, we present some po¬ 
larisation simulations for aspherical but axially symmetric configu¬ 
rations. Our results for a prolate morphology with a pole-to-equator 
density ratio of ~ 3 yield a maximum < 2% continuu m polarisa¬ 
tion c ompatible with the value observed for SN 2010jl jPatat et al.l 
l201lh . In this context, the asymmetry may stem only from the 
CSM, although both inner and outer shells may be asymmetric. In 
future work, we will extend the HERACLES simulations to 2-D to 
study large scale asymmetries and to 3-D to study the stability of 
the dense shell. 
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Figure 21. Velocity profiles from the HERACLES simulations for models 
Xel0m6 (top), Xe3m6r (middle), and Xm6 (bottom). The epochs shown 
con-espond to 0.01,1, 2, 3, 4, 5, 6, 7, 8, 10, 20, 30, 50, 100, 150, 200, 250, 
and 300 d after the onset of the interaction. Notice the pre-acceleration that 
takes place ahead of the shock. 
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Figure 22. Left: Maximum light spectra for the interaction models with different ejecta kinetic energy and/or CSM mass/density (see Table^for additional 
details). Note the strong shift in colour for interaction models with higher ejecta kinetic energy. Right: Same as left, but showing the H« line region, in velocity 
space. Interestingly, even for very large CSM densities (and thus mass loss rates), the maximum extent of the electron-scattering wings is 2000-3000 km s“^. 
The much broader wings of Ho in SN 2010jl as time progresses suggest that they do not stem exclusively from scattering with electrons but must instead 
contain a contribution from Doppler-shifted line photons (i.e., emitted from the fast moving material in the interaction region). 
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